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Abstract 


We design a fast technique for fitting cubic Bézier curves to the boundary 
of 2D shapes. The technique is implemented by means of the Nelder—Mead 
simplex procedure to optimize the control points. The natural attributes 
of the Bézier curve are utilized to discover the initial vertex points of the 
Nelder—Mead procedure. The proposed technique is faster than traditional 
methods and helps to obtain a better fit with a desirable precision. The 
comparative analysis of our results describes that the introduced approach 
has a high compression ratio and a low fitting error. 
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1 Introduction 


Capturing the outlines of 2D objects is a substantial topic in computer aided 
geometric designs (CAGD), computer graphics, as well as vision and imaging; 
see [1, 2,6, 7, 9-12, 15, 20-25, 32, 38, 39]. Curve fitting with Bézier curves is a 
traditional problem in this process. The goal of a curve fitting method is to 
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detect a collection of the control points that can precisely indicate the given 
target shape. 

In the literature, many methods have been suggested to represent out- 
line objects [3, 11, 21,29, 34]. A vast majority of these methods pertain to 
an optimization problem that finds an appropriate curve for the data gen- 
erated by the outlines of 2D shapes. Most common algorithms for solving 
this model are not suitable optimization algorithms. Since they involve com- 
plex computations, they cannot be used for real time applications. Itoh and 
Ohno [11] exerted the least square fitting to approximate cubic Bézier curves 
without fixing the end points of the curves. Plass and Stone [21] proposed 
an iterative procedure for fitting a parametric piecewise cubic polynomial 
curve with an selective endpoint and tangent vector specifications. Srafraz 
and Khan [28] employed the least square method to certify an appropriate fit. 
Their method contains several phases, such as extraction of boundaries, ex- 
ploration corner points and break points, and fitting curves. In another study, 
the same authors added reparameterization steps to ameliorate the efficiency 
of the fit [29]. Sarfraz and Razzak [34] employed a generalized Hermite cubic 
spline to capture the outline of digital character images using characteristic 
points and the least squares method. The enhanced Bézier curve scheme pro- 
posed by Sohel et al. [36] reduces the distance among the Bézier curve and 
its control polygon without extra computational complexity. An object cod- 
ing technique using Bézier curves was introduced by Masood and Haq [16]. 
They determined the control points by searching along the endpoint tangents. 
Sarfraz and Masood [30] determined the appropriate position of the control 
points by using the natural attributes of cubic Bézier curves. Masood and 
Sarfraz [18] presented an outline capturing scheme using Bézier cubic approx- 
imation. Their method operates in two phases to find intermediate control 
points. In Phase 1, the location of the detected control points gets closer 
to that of the original control points. Phase 2 is exerted to hit the target 
location very exactly. Masood and Sarfraz [17] used the properties of cubic 
Bézier curves and analysed the control points spread (CPspread) to subdi- 
vide complex segments into two or more segments. Some other techniques 
based on the search algorithm include evolutionary algorithm [35], simulated 
annealing approach [26,27], genetic algorithm [32,33], and wavelets [37]. 

In the above mentioned methods, the control points are specified by the 
least square algorithm or using the properties of cubic Bézier curves. These 
processes are computationally expensive and not appropriate for practical 
applications. Our contribution in this study is to introduce a scheme to avoid 
these time consuming operations and to guarantee a desirable equivalence 
between compression ratios and fitting errors. We perform the Nelder-Mead 
algorithm to find the intermediate control points of the cubic Bézier curve. 
The initial situation of the control points is specified using extracting some 
effective virtues of the Bézier curve. 

The rest of the paper is arranged as follows. Section 2 reviews the seg- 
mentation of an outline by corner detection. Section 3 introduces the Nelder— 
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Mead algorithm. The process of determining the control points for the cubic 
Bézier curve approximation using the Nelder-Mead simplex method is de- 
scribed in section 4. Section 5 illustrates the experimental observations and 
compares them with those of other methods. Finally, the paper is closed up 
with a conclusion in Section 6. 


2 Outline segmentation 


Outline segmentation is used to divide the shape outline into small part and 
simplify the curve fitting procedure. We use the corner points to partition 
the outline into disparate segments from the natural break points. Authors 
have suggested various corner detection algorithms in the literatures [3,4,31]. 
These algorithms do curvature analysis with numerical techniques. In this 
paper, the algorithm designed by Sarfraz et al. [31] is used for the corner 
detection, because this algorithm is accurate, effective, and robust to noise. 
The method is concisely explained in this section (readers are referred to [31] 
for details). The algorithm detects corner points in two steps. Candidate 
corner points are discovered from the outline data points in the first step. If 
all the boundary points are Q;,1 <i<n, then the boundary point Q; can 
be given as 


if (¢@+L)<n, 
then Q, = Qitz, 
else Qk = Qi41)-n; 


where L represents the length parameter and preserves of the shape scaling 
and resolution. The default value of L is 14. The perpendicular distance d; 
from point Q,(z, y) to the direct line joining the points Q; (x,y) and Q(z, y), 
can be given as follows: 


if mz = 0, 
then dj = |Qj.0 = Qi,z|; 
lQjy — MQj,0 +MQie — Qiyl 


else dj = Jeet 
m Qky — Qi 

where m = —4 = =" ag 
Mea Qk,x _ Qin 


The point @,; is specified as a volunteer corner point if its perpendicular 
distance (d;) from the direct line Q;Q, goes beyond D. The local sharpness 
and the opening angle of the corners are investigated by the distance param- 
eter D. It also checks the incorrect determination of corner points that may 
occur due to noise and disarray. The default value of D is set to be 2.6. The 
new direct line through increasing both 7 and & detects the next volunteer 
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corner points. In the second step, the superfluous corner points are deter- 
mined. A superfluous corner point is one of that any other volunteer with a 
higher value of d; is in the domain R. In the other words, a candidate corner 
point will stay if it has the greatest value of d; between the R number of 
points on its both sides. The default value of R is equal to length parameter. 
In the present study, we utilize this method at its default values for corner 
detection. All steps of this corner detection method is shown in Algorithm 
lg 


Algorithm 1 The corner detection algorithm 


Set values of L, D, and R; 
Detect boundary points {Q;}"_, of 2D object using the Canny edge de- 
tector; 
First step: detect candidate corner points; 
for i=1tondo 
if (¢+L)<nthen 


Qk = Qitz, 
else 

Qk = Qiitb)—n3 
end if 


for 7 =itokdo 
Calculate d; the perpendicular distance from point Q; to the 
straight line joining points Q; and Qx; 
if dj > D then 
Add Q, to candidate corner point; 
end if 
end for 
end for 
Second step: detect superfluous corner points; 
A candidate corner point will stay if it has the greatest value of d; between 
the R number of points on its both sides; 


3 The Nelder—Mead simplex method (NMSM) 


The Nelder—Mead algorithm is a numerical optimization procedure for solving 
unconstrained optimization problems in a multidimensional space; see [5, 14, 
19]. This method is one of the enormously famous derivative-free techniques. 
It is a direct search algorithm that only needs a numerical evaluation of the 
objective function at a finite number of points per iteration. This method is 
planed for prevalent unconstrained minimization problems, such as nonlinear 
least squares and nonlinear simultaneous equations; see [8, 13]. 
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The problem discussed in this study is the following unconstrained opti- 
mization problem 


min f(z) «€ER”, 


where f is a nonlinear function from R” into R and xz € R”. For the function 
f(a), the Nelder-Mead method begins with n+ 1 vertices as yi, ya,---;Yn+1 
and then evaluates the objective function value of each vertices. We assign 
to yi as the best vertex where f(x) is the lowest and to y,+41 as the worst 
point where f(x) is the highest. The most common Nelder—Mead iterations 
perform a succession of primary geometric transformations including reflec- 
tion, expansion, contraction, and shrinkage to find a better point and make 
it replace the worst point. Geometric transformations are controlled by four 
parameters including coefficients of reflection p, expansion y, contraction 7, 
and shrinkage o. The introduced parameters should assure the following 
consternations: 


p>0, x>1, 0<y<1, and 0<o<1l. 


The common values, exerted in this paper, are 


1 1 
p=1, x=2, ¥ rt and oS 


We set ¢ = 10~!° for the termination criterion in the Nelder-Mead method. 
The principal foundation of the Nelder-Mead method is described in Algo- 
rithm 2. Readers are referred to [8,19] for a detailed study of the Nelder-Mead 
method. 


4 Curve approximation 


We partition the outline of a 2D shape into curve segments based on its 
corner points. It means that all the outline points between two consecutive 
corner points constitute one curve segment. In this literature, cubic Bézier 
curves have been used to fit each segment separately using the Nelder-Mead 
algorithm. 


4.1 Bézier curve 


The Bézier curve is a parametric curve C(t) that widely used to describe and 
model curves and surfaces in computer aided designs and computer graphics. 
The Bézier curve of degree n is given as 
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Algorithm 2 The Nelder—Mead algorithm [8] 


1. Choose an initial n + 1 vertex points {y1, y2,---,Yn+i} and choose a 
stopping criterion €; 
2. Order. Order and re-label the n + 1 vertices to assure f(y) < 


f(y2) S++ S FYn41)3 
3. Reflect. calculate the reflection point y, by yr = 9 + p(¥ — Yn+1); 


where y = yo is the centroid of the n best points; 
n 


if f(y) < fr) < fn) then 
replace y,+41 with the reflected point y, and go to Step 7; 
end if 
4. Expand. 
if f(y-) < f(y) then 
calculate the expansion point ye by ye = 9 + x(Yr — 9); 
end if 
if f(ye) < f(yr) then 
exchange Yyn+1 with ye and go to Step 7; 
else 
exchange Yn+1 with y, and go to Step 7; 
end if 
5. Contract. 
if f(yr) = f(yn) then 
a contraction is performed among y and the better of yn+1 and y;; 
end if 
a. Outside. 
if f(Yn) < f(Yr) < f(Yn41) then 
apply an outside contraction: compute Yoo = ¥ + ¥(Yr — 9); 
end if 
if f(Yoc) < f(yr) then 
exchange Yn41 With yo- and go to Step 7; 
else 
go to Step 6 (perform a shrink); 
end if 
b. Inside. 
if f(yr) = f(Yn+i) then 
apply an inside contraction: compute yjc = ¥ + Y(Yn41 — 9); 
end if 
if f(Yic) = f(Yn+i) then 
exchange Yn+41 with y;- and go to Step 7; 
else 
go to Step 6 (perform a shrink); 
end if 
6. Shrink. Measure the n new vertices 
y =yit+o(yi — yi),t =2,...,n+1. 
exchange the vertices y2,...,Yn+1 with the new vertices y,..., Yh413 
7. Termination Criterion. Arrange and re-label the vertices of the 
new simplex as yj, Y2,---,Ym+1 such that f(yi) < f(y2) <--> < f(Yn41); 
if f(Yn+i) — f(y) < € then 
stop; 
else 
go to Step 3. 
end if 
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Ci) =>" PB Bint), 0<¢<1, 
1=0 


where P; refers to the control points and B;,,(t) is the Bernstein basis poly- 
nomials of index 7 and degree n expressed as 


n 


Bin(t) = ( Jea-or, = D,e say: 


t 


From the defining equation of the Bézier curve, it can be seen that the prop- 
erties of Bernstein polynomials are passed on to the Bézier curve; see [20]. 
These properties are as follows: 


e Partition of unity. 

e Affine invariance. 

e Convex hull property. 

e Endpoint interpolation. 


e Value of Bernstein polynomials, for all 0 < ¢ < 1, does not appertain 
to the location of any control point(s); see [17]. 


e If the location of any control point is specified, its efficacy all over the 
curve can be deprived; see [17]. 


4.2 Curve fitting with the cubic Bézier curve 


Let {X;}_, C R? denote an ordered set of contour points of a segment. Our 
goal is to find control points P; (i = 0,1, 2,3) so that the cubic Bézier curve 


C3(t) = Po Bo,3(t) + P, By 3(t) + P2Bz,3(t) + P3B3.3(t), 0 <t< 1, 


can be a good representation of {X;}+_,. For this purpose, we minimize the 
linear least squares fitting error, E’, specified as the sum of the squares of the 
deviations: 


N N 3 2 
E =o |lCs(tu) — Xell? = >|] 0 PBi alte) — Xe 
k=1 k=1 i=0 


where the parameter values t; are assigned to X; using the uniform parame- 
terization method; see [20]. The control points Py and P3 are the two corner 
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points of the segment. So, the effect of Py and P3 is removed from a given 
cubic Bézier curve as 


N 
E =o |ICy(te) — X4IP, (1) 
k=1 
where 
C3(t) = Pi Bi,3(t) + P2Ba,3(t) 
and 


Xj, = Xx — (PoBo,3(te) + P3Bs,3(te)), &=1,2,...,N. 


Therefore, we have to find the control points P, and P2 by minimizing the 
problem defined in (1). 

Let P; = (P2,, Py,),i = 1,2, and write all the intermediate control points 
in one vector: 


P, 
Py, 
P,, 


P. 


y2 


1 


With this vector, the fitting problem is formulated as 
N 
min) ||C4(éx) — X41. (2) 
k=1 


That is, there are four parameters to estimate. We employ the Nelder-Mead 
method to solve the least square fitting problem defined in (2) and find these 
parameters. 


4.3 Initialization 


Before Algorithm 2 is implemented, a starting point should be selected for 
the intermediate control points. The initial vertex points can be chosen 
arbitrarily, but using the appropriate initial situation of the control points 
increases the speed of convergence and quality of solutions. Duo to x € R# 
in the objective function (2), five vertices have to be chosen. The details on 
the selection of the initial vertex points are explained as follows. 

Let 


P, By,3(t) + P2Bo,3(t) = C3(t), (3) 


Capturing outlines of generic shapes with cubic Bézier curves ... 11 
where 
C3(t) = C3(t) — PoBo,a(t) — P3Bs,a(¢). 
Based on the properties of the cubic Bézier curves, we get 
B = B,3(0.5) = Bo 3(0.5). (4) 
From equations (3) and (4), we have 
Pi + Pp =Cs, (5) 


where 


~ _ C3(0.5 
G, = 305), 


By solving equations (3) and (5), we have 


R= C3(t) — Bo3(t)C3 
* “Bi 3(t) — Bas(t) ’ 


Py = C3 — Py. 


Suppose P,' and P,' are the computed intermediate control points at 
(0 SUS 1, tj # 0, 0.5, L,), and 


Pig C3(ti) — Bo,3(ti)C3 
1 = ) 
By 3(ti) — Bo,3(ti) 


Pp, =C;-— P,*. 
where 7 = 1,2,...,5. 
Let Pi = (Pi,, Pi) and Pj = (P3,, Pi,). The vector x; can be identified 


w) @2? 
as 
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a= ||, 1=1,2,...,5. 


Hence, these vectors can be used for the initial vertex points of the Nelder— 
Mead method. 

In order to illustrate the curve fitting using the Nelder-Mead algorithm, 
we consider a curve as shown in Figure l(a). Figure 1(b) shows the fitting 
of a given curve (black) using a cubic Bézier curve (red). The initial vertex 
points are marked by e, and the intermediate control points computed by the 
Nelder—Mead algorithm are marked by «. The fitting error versus the number 
of iterations for this curve with different initial vertex points is demonstrated 
in Figure 2. As shown in Figure 2, the Nelder-Mead method based on the 
proposed initial vertex points converges faster than the one based on random 
initial vertex points. 


- 7 


(a) Given curve. (b) Cubic Bézier approximation. 


Figure 1: Fitting a given curve (black) using a cubic Bézier curve (red) with 
the proposed algorithm (NMSM). 


4.4 Segment subdivision 


When curve fitting is done by means of the proposed method (i.e. NMSM), 
there is a possibility for a generated complex curve not be acceptable and 
for the deviation error to be high. In such a case, the segment subdivision is 
applied to decrease the fitting deviation. In this paper, the complex segment 
of the outline is partitioned into two parts at the point of the maximum 
deviation error if the maximum deviation error oversteps the given error 
threshold limit (ETL). The process of subdivision will continue constantly 
and recursively until the fitting error reaches a value under the given error 
threshold. An example of subdivision into five curves is demonstrated in 
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Proposed initial vertex points 


——Random initial vertex points 


5 10 1520 25 30 3540 45 60 
iteration 


Figure 2: Error vs iteration. 


Figure 3. Figure 3(a) is the curve fitting without segment subdivision and 
Figure 3(b) is curve fitting after segment subdivision. The black and red 
lines display the original and the approximated curves, respectively. The 
subdivided points are marked by M@. All the steps of the outline capturing 
system are shown in Algorithm 3. 


(a) Without subdivision. (b) After subdivision. 


Figure 3: Curve fitting with recursive subdivision. 
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Algorithm 3 The outline capturing system 


Get a digitized image and choose an error threshold limit (ETL) ; 
Extract an outline; 
Detect the corner points and divide the outline into segments; 
for each segment do 
Compute the initial vertex points; 
Perform curve fitting over the segments using the Nelder—Mead algo- 
rithm; 
Compute the maximum deviation error (MDE); 
if MDE < ETL then 
break 
else 
Find the subdivision point and partition into two segments; 
end if 
end for 


5 Experimental results and comparative study 


The proposed algorithm explained in the Section 4 has been used to some 
explanatory examples corresponding to generic shapes. Then, the outputs 
have been compared with those of algorithms [17, 29,34]. In order to allow 
a fair comparison, we have tested the same shapes of the criterion as in 
other papers. In this study, the results are measured based on the following 
parameters: 


e Number of segments: The outline of the shape is partitioned into a 
number of segments using a corner detection algorithm and subdivision 
points. An increase in the number of segments, leads to an increase in 
the number of control points. Therefore, the number of segments must 
be decreased. 


e Compression ratio (CR): This is a significant criterion in the evaluation 
of the amount of compression accomplished by the method. A large 
compression ratio is favorable. This criterion indicates the ratio of the 
number of data points in an actual outline (n) to the obtained data 
points in curve fitting (nDP). It can be computed as: 
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e Average error: This parameter refers to all the errors generated in the 
fitted outline of the shape and is given as: 


n 
1 
Average error = — y C4, 
n 
i=1 


where e; is the Euclidean distance between the ith point of the outline 
and the corresponding point of the parametric curve. 


e Maximum deviation: This parameter calculates the maximum devia- 
tion of the fitted outline from the original shape. It is calculated by the 
following equation: 


. . . n 
Maximum deviation = max{e;} 
i=1 


e Computation time: This parameter indicates the amount of time taken 
by the algorithm for capturing the outline of the shape. It appertains 
to the performance approach and the processor used. 


The algorithm has been performed and compared for Figures 4 and 5, 
which are Arabic words “Sabr”, and “Kanji” characters, respectively. Figure 
4(a) shows the image of the Arabic word “Sabr” and its boundary is given 
in Figure 4(b). Figure 4(c) shows the boundary along with the corner points 
marked with a square (Ml). Figure 4(d) shows the captured outline based 
on our approach, and the subdivision points are denoted with a circle (e). 
The outline computed by means of the method presented by Masood and 
Sarfraz [17] is shown in Figure 4(e). Figures 4(f) and 4(g), respectively, 
represent the outline evaluated by Sarfraz and Razzak [34] at a threshold of 
3 and 2. The quantitative comparison of the different methods in Figure 4 
can be observed in Table 1. 

The results obtained for the “Kanji” character are presented in Figure 5 
and tabulated in Table 2. The image of “Kanji” character and its boundary 
are shown in Figures 5(a) and 5(b), respectively. Figure 5(c) shows the 
boundary along with the corner points which are marked with a square (MH). 
The results of the presented method for this shape are given in Figure 5(d). 
The outline captured with algorithm [17] is shown in Figure 5(e). The result 
obtained by Sarfraz and Khan [29] can be seen in Figure 5(f). The outline 
generated by means of method [34] is given in Figure 5(g). It should be noted 
that the results of Tables 1 and 2 are achieved by different systems. 

Through a visual inspection of the results in Tables 1 and 2, one may 
come to the following conclusions: 
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Table 1: Quantitative comparison for the Arabic word “Sabr”. 

Algorithm Number of | Compression | Maximum Average Computation 
segments ratio deviation error time 
Proposed (NMSM) 28 46.43 1.92 0.44 0.79 
Masood and Sarfraz [17] 27 48.15 2.21 1.38 0.83 
Sarfraz and Razzak [34] at 7 = 3 66 19.70 1.71 1.07 3.79 
Sarfraz and Razzak [34] at 7 = 2 98 13.27 1.39 0.91 2.58 

Table 2: Quantitative comparison for the “Kanji” character. 

Algorithm Number of Compression Maximum Average Computation 
segments ratio deviation error time 
Proposed (NMSM) 33 64.81 1.53 0.40 0.55 
Masood and Sarfraz [17] 33 64.81 1.55 0.73 0.59 
Sarfraz and Khan [29] 31 69 3.78 1.92 3.61 
Sarfraz and Razzak [34] 66 33.42 1.40 1.13 3.83 


The method proposed in this study protects an appropriate equilibrium 
among the compression ratio and the fitting error. 


The proposed method avoids time consuming procedures, such as least 
square fitting, and uses one of the popular derivative-free free operations 
(ie., the NelderMead simplex method ) that is simple and fast. 


The average error, maximum deviation, and the computation time are 
lesser than the results obtained by other methods in the current liter- 
ature [17, 29, 34]. 


The number of segments generated by method [34] is very high because 
the corner detection and the subdivision algorithms are suboptimal. 
The suboptimal corner points are indicated by arrows in Figures 4(f), 
4(g), and 5(g). This method also suffers from a long computation time. 
This is due to the use of least squares fitting and control parameters. 


Computation of the outline by means of method [29] takes too long 
and has errors. The least squares curve fitting and the noise filtering 
procedure are the main causes for these disadvantages. It is shown with 
an arrow in Figure 5(f). 


The average error with algorithm [17] is higher than that in the pro- 
posed method, because this method does not use the standard opti- 
mization process. 
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(d) Proposed method (NMSM). (e) Masood and Sarfraz [17]. 
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(f) Sarfraz and Razzak [34] atr = 3.  (g) Sarfraz and Razzak [34] at 7 = 2. 


Figure 4: Captured shape of the Arabic word ”Sabr” with different methods. 
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(d) Proposed method (NMSM). 


(f) Sarfraz and Khan [29]. 
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(b) Outline of original (c) Boundary with cor- 
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(e) Masood and Sarfraz [17]. 
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(g) Sarfraz and Razzak [34]. 


Figure 5: Captured shape of the Kanji character with different methods. 


6 Conclusion 


In this study, we introduced a novel outline capturing scheme for 2D shapes 
based on the Nelder—Mead simplex method. The Nelder—Mead simplex algo- 
rithm is straight to perform, easy to calculate, and does not call for much time 
to compute the optimization problem. The unique feature of this method is 
that it does not require to consider time consuming least squares fitting. The 
initial situation of the control points is obtained through the properties of 
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the cubic Bézier curve. It increases the speed of convergence and the quality 
of curve fitting. As a result, the introduced approximation scheme is much 
faster than traditional approaches and helps to capture a better fit with a 
desirable precision. Through a comparison of our method with previous ap- 
proaches and considering the simulation results, it emerges that the method 
has such advantages as appropriate equivalency between the compression ra- 
tio and the fitting error, low deviation error, and low computation time. 
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